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RESUMO 


O fenômeno de autoexcitação em máquinas síncronas tem por razão os efeitos de ressonância 
elétrica aos quais estes equipamentos ficam sujeitos em função da presença de cargas indutivas 
associadas a cargas capacitivas. O resultado deste fenômeno pode ser acompanhado de fortes 
sobretensões tanto no campo quanto no estator da máquina. Esse fenômeno também pode ser 
agravado durante as oscilações no sistema interligado uma vez que o processo de ressonância 
tem a ver com a frequência elétrica envolvida. Este artigo tem por objetivo apresentar uma análise 
sucinta do ponto de vista analítico para, posteriormente, calcular a faixa de risco operativo 
considerando o fenômeno da autoexcitação. Serão utilizadas as chamadas transformações 
ortogonais de Park, também conhecidas como transformações de Blondel, seu primeiro 
idealizador. Na realidade a transformação a ser utilizada será uma variante daquela que foi 
sugerida por Anderson [1]. 
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1. Explanação Geral Acerca De Máquinas Síncronas 


Ao se iniciar o estudo relativo à máquina síncrona, pode-se apresentar de início uma comparação 
direta entre as máquinas síncronas e outras máquinas rotativas. 
Com relação às máquinas assíncronas ou máquinas de indução, outro tipo físico de máquina elétrica, 


pode-se afirmar que: 


a) as máquinas síncronas, quando operam como gerador, fornecem, ao mesmo tempo, potência ativa 
e potência reativa, já o gerador assíncrono fornece apenas potência ativa, continuando como no 


motor assíncrono, a consumir potência reativa; 


b) um gerador assíncrono consome em média, por característica, uma potência externa de 
aproximadamente 40% de sua potência nominal no processo de excitação o que inviabiliza sua 
utilização já que o processo de excitação de um gerador síncrono exige apenas uma potência da 


ordem de 1% para menos de sua potência nominal. 


A máquina síncrona pode ser representada, esquematicamente, pela figura 1. 
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Figura 1 — Esquema Básico de uma Máquina Síncrona com um Par de Pólos e Ligação do Estator 
em Estrela. 


Pode-se notar pela figura 1 que existem dois eixos, um na direção do pólo e outro em quadratura 
com este. O eixo na direção do pólo é chamado de eixo direto, ou apenas d, enquanto o outro é 
chamado de eixo em quadratura ou apenas q. Apenas lembrando que nesse caso, por questões 
didáticas, o estator se refere à parte fixa da máquina e o rotor à parte que se move com velocidade 
w em relação ao estator. 

Observar que este sistema de eixos d-q, idealizado, gira na mesma rotação do rotor, ou seja, com 
velocidade q. 

Como o torque mecânico é transferido através da rotação do rotor via eixo de acoplamento rotor/fonte 
de energia primária, quando em gerador síncrono, ou via eixo de acoplamento rotor/carga mecânica, 
quando em motor síncrono, existe um sincronismo entre o campo magnético de excitação e o torque 
mecânico que é desenvolvido neste eixo de acoplamento. Isso quer dizer que ambos os agentes, 
campo magnético e torque mecânico, deverão girar na mesma velocidade, portanto sincronizados 
em regime permanente. Esta é a razão de a máquina se chamar máquina síncrona. 

Nota-se, também, que o sistema de eixos d-qg gira sincronizado na frequência de rotação do conjunto 
girante uma vez que se trata de eixos fictícios idealizados sobre o rotor da máquina. A idealização 
desse sistema de eixos vem ligada às facilidades de análise quando todas as grandezas, mesmo 
aquelas localizadas no estator, são referidas ao rotor, portanto aos eixos d-g. Estas transformações 
são clássicas no estudo da máquina síncrona e são conhecidas como transformações de Blondel, 
primeiro estudioso que as idealizou. Logo após, Park, Doherty, Nickle, entre outros estudiosos, 
apresentaram trabalhos nesta área. Nesse livro essas transformações serão chamadas apenas de 
transformações d-g-0. 

O estudo da máquina síncrona, considerando-se as transformações de Blondel, será então 
desenvolvido de maneira fundamentada uma vez que sem esse conceito é difícil compreender 


perfeitamente o processo de regulação de tensão da máquina síncrona. 


1.1. Enrolamento de Campo 


O enrolamento de campo que envolve cada pólo cumpre o papel de, a partir de uma fonte de corrente 
contínua, que pode ser externa ou auto-alimentada, produzir o campo magnético indutor responsável 
pela excitação da máquina síncrona. 

Em uma máquina síncrona com um par de pólos, por exemplo, nota-se que o fluxo magnético gerado 
em um dos lados do eixo polar é norte enquanto que do lado oposto o fluxo magnético gerado é sul. 
Em outras palavras, trata-se de um grande eletrofmã. O enrolamento em si é constituído por espiras 
na forma de barras de cobre eletrolítico numa disposição retangular que envolve cada pólo de modo 


concentrado. 


1.2. Estator 


O estator da máquina síncrona, na sua forma padrão, constitui a parte estática que é varrida pelo 
fluxo magnético contínuo gerado no rotor. É formado pela carcaça, pelo núcleo magnético e pelos 
enrolamentos de fase da máquina. 

O núcleo do estator, semelhantemente ao núcleo do rotor, é confeccionado de chapas de aço 
silicioso, ou seja, uma liga de aço com silício que apresenta alta permeabilidade magnética e baixas 
perdas, isoladas umas das outras e montadas segundo um alinhamento definido previamente através 
de pinos guia e presas por tirantes. O conjunto é cnumbado na estrutura civil da instalação elétrica. 
O enrolamento do estator é constituído de bobinas cujos condutores encontram-se alojados em 
ranhuras na sua parte interna. A estrutura final do anel magnético do estator termina em uma 
armação circular tendo na sua parte externa a carcaça e os trocadores de calor e em sua parte interna 
uma superfície dentada com as ranhuras que comportarão as bobinas. O estator da máquina 


síncrona também é conhecido como armadura da máquina. 


1.3. Sistema de Excitação 


Este sistema é responsável pela produção do fluxo magnético indutor a partir de uma fonte de 
corrente contínua podendo ser externa ou auto-alimentada que, no caso de grandes máquinas, se 
constitui de um controlador automatizado e que em pequenas máquinas pode ser composto de um 
sistema reostático controlado manualmente. Nas grandes instalações, esse sistema torna-se 
complexo e cumpre importantes ações que se refletem na própria estabilidade do sistema elétrico de 
potência, o SEP. 

A tensão de corrente contínua de alimentação é levada ao enrolamento de campo através de escovas 
que se encontram em contato permanente com os anéis coletores. 

Além daquilo que foi descrito até aqui, é importante lembrar que as máquinas síncronas se diferem 
quanto à construção do rotor e, por conseguinte, do enrolamento de campo. 


2. Equações Da Máquina Síncrona 


2.1. Geral 


A máquina síncrona caracteriza o elemento representativo do processo sobre o qual vão atuar os 
órgãos de comando e controle para a obtenção do produto final, qual seja, a potência elétrica 
despachada a partir da potência motriz de entrada fornecida pela turbina. 

Para o equacionamento da parte relativa à regulação de tensão, faz-se necessária a utilização da 
representação da máquina síncrona através de um equacionamento tal que venha a facilitar as 
análises e os estudos envolvidos. 

O procedimento normal é a aplicação das chamadas transformações d-g-0 que busca a conversão 
de grandezas estatóricas em grandezas rotóricas através da transformação de grandezas gerais em 
grandezas pseudo-estacionárias referidas a um sistema de eixos fixado sobre o rotor e que gira na 
mesma velocidade deste, velocidade síncrona, fazendo com isto que estas grandezas independam 
da rotação. 

O fato consiste na criação de três eixos chamados d-g-0, ou mais propriamente, eixo direto, eixo em 
quadratura e eixo zero, sendo que os eixos d e q perfazem entre si um ângulo de 90 graus elétricos 
e se movimentando na velocidade síncrona, sendo o eixo zero um eixo estacionário simplesmente 
utilizado nas representações em que as condições desequilibradas são verificadas. 

O processo visa à projeção das grandezas dos eixos estatóricos, aqui designados por a-b-c, segundo 
uma matriz de transformação particular, sobre aqueles eixos d-g-0, do rotor, e a partir daí obter o 
equacionamento de toda a dinâmica da máquina síncrona. 

Com o estabelecimento do equacionamento básico relativo às indutâncias, segundo as influências 
próprias e mútuas dos enrolamentos envolvidos, serão equacionadas as funções de transferência 
das partes implicadas buscando-se com isto a definição do modelo da máquina síncrona. 

Em cima da modelagem e da análise, detalhadas previamente, será possível a realização de 
diagnósticos sobre o comportamento operativo da máquina síncrona sob determinadas condições. 


2.2. À Transformação d-q-0 


Como já foi citada anteriormente, a idéia da utilização da transformação de grandezas estacionárias 
em grandezas pseudo-estacionárias teve seu embasamento matemático formulado na utilização de 
uma matriz de transformação, um método específico de tratamento matemático para essa questão. 

A transformação a ser aqui utilizada se trata da transformação específica proposta por W A Lewis, 


ver referência 8. 


Primeiramente será representada uma máquina síncrona composta de um par de pólos conforme 
está mostrado na figura 1. 

Os eixos d e q giram com a velocidade síncrona do conjunto no sentido anti-horário. O eixo d 
encontra-se em fase com a direção do pólo enquanto o eixo q encontra-se em fase na direção 90 
graus elétricos atrasados do eixo d, considerando-se o sentido de giro. 

O diagrama fasorial para a máquina síncrona, anteriormente mostrado na figura 1, é apresentado na 


figura 2 considerando-se apenas a fase a do estator. 
A equação básica da transformação d-g-0 é a seguinte: 


Gago = P Gabe (001) 


Onde: 


G: matriz das grandezas no sistema indicado; 
P: matriz da transformação. 


Sentido de Rotação 


s— () 


Lp Eixo d 


Figura 2 — Diagrama Vetorial da Máquina Síncrona para a Fase a 


Na figura 2 considerar que o ângulo elétrico de giro do sistema seja dado pela seguinte expressão: 


B(1)=S+wt+T (002) 


onde w é a rotação ou frequência do conjunto girante, geralmente tomada como um valor fixo, 
síncrono e sob controle automático em função da presença do regulador de velocidade. 

Em relação à figura 2, deve-se notar que essa representação expressa a posição momentânea dos 
vetores de tensão da máquina síncrona, mormente o vetor de tensão terminal da fase a, vetor va em 
relação ao vetor da tensão de excitação Ea que se encontra sobre o eixo q. Deste modo, num dado 
momento em que a localização do rotor fica definida, o vetor da tensão terminal relativa à fase a está 
localizado exatamente a um ângulo 6 radianos elétricos atrás do eixo d' cuja flecha positiva indica o 
sentido da fmm perpendicular ao plano de uma espira hipotética que concentraria toda a corrente 
gerada na fase a pelo processo de conversão eletrodinâmica da máquina elétrica. A equação 002 é 
bem clara nessa interpretação. 


Nota: A partir desse ponto, no sentido de dar mais clareza e simplicidade ao equacionamento que 
será um tanto quanto trabalhoso, nas designações de grandezas variantes no domínio temporal 
será suprimida a denotação (t) como indicação de dependência da variável tempo. 

Assim, a equação 002, de modo mais simplificado, tomará a seguinte forma: 


g-orwt+T 
2 


Considerando- se a matriz P, da transformação desejada, vem que: 


1 1 1 

Vo Ha Vi 

p= fé cos6 cos (60-21) cos (6+2rm) 
3 3 3 

2r) 
3 


sin9 sin(o-Sm) sin (0 + 


(003) 


A transformação P, aqui utilizada, trata-se de uma transformação ortogonal, ou seja, aquela que 


obedece à seguinte operação quando a inversa da matriz original é exatamente igual a sua 


transposta: 
p!=p” (004) 
Portanto, para as correntes, para as tensões e para os fluxos magnéticos enlaçados, grandezas 


variantes no tempo e associadas à armadura da máquina síncrona e refletidas no sistema d-g-0, 


pode-se escrever, de forma matricial resumida, que: 


idgo =P labe (005) 


Vdgo = P Vabe (006) 

Adgo = P Aabc (007) 
A matriz inversa da matriz P de transformação é a seguinte: 

par cos6 sing 


—1 2 1 2 e ne 
P É Ed cos (8 Guê, sin(O ué, 


Ze cosp+ám) sing + Em) (008) 


A transformação escolhida garante a invariância da potência em ambos os sistemas a-b-ce d-g-0 e 
para comprovar isso bastará considerar a expressão da potência instantânea, no domínio do tempo, 
deliberada por uma máquina síncrona operando como gerador. 

Assim, vem que, levando-se em conta as tensões e correntes instantâneas de cada uma das fases 


do estator da máquina síncrona, variantes no tempo: 

pP=Vaia+Vbib+Veic (009) 
Onde:p: potência instantânea total gerada; 
Va bc: tensões instantâneas de estator por fase; 


lab: Correntes instantâneas de estator por fase. 


De outra maneira, a equação 009 pode ser reescrita numa forma mais adequada para as operações 


matriciais futuras. Deste modo, tem-se que: 


p=[vabe]! Label (010) 


Notar que as tensões e as correntes instantâneas de estator da máquina síncrona se distribuem na 
forma de uma matriz coluna, assim como tais entidades são denominadas na teoria. Isso quer dizer 


que: 


na (011) 


[vabe] =|Vb |=Vabc 
Vc 


ia 
Liabel= ib = iabe (012) 


lc 


Notar que na equação 010, a matriz coluna das tensões de fase teve de ser transposta para que a 
operação de multiplicação matricial tivesse consistência. 
A partir das equações 005 e 006, podem-se escrever as seguintes equações, respectivamente para 


as correntes e as tensões de fase da máquina síncrona: 
np; e o. 
labe = Idgo Vabce =P Vdgo 


Levando-se as expressões anteriores, respectivamente para as correntes e para as tensões, na 


equação 010, vem que: 
p= La "vago te 'iago |- [ano] [” É Es lidgo (013) 


Uma importante observação relativa à transformação via matriz P é que ela se trata de uma 
transformação ortogonal, valendo então a relação da equação 013 anterior e as observações 
concernentes a esta. 

A partir da matriz de transformação P, tem-se, para a equação 013 referente à potência total 


instantânea transferida pela máquina síncrona, que: 
li 2 


Através do cálculo matricial anterior fica comprovada a invariância da potência na transformação a- 
b-c — d-g-0. Essa propriedade será fundamental para o estudo da transformação Gago - ? Gabc e vice 
versa. 

Uma observação interessante: quando se trabalha com sistemas trifásicos equilibrados, as projeções 
das grandezas sobre o eixo zero tornam-se nulas uma vez que a existência desse eixo tem a ver 


apenas com a questão de um possível desequilíbrio entre as grandezas trifásicas. 
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2.3. Equação Geral dos Fluxos Enlaçados 


Na figura 1, embora muito simples, podem-se verificar os enrolamentos que interferem no processo 
de trabalho de uma máquina síncrona. 

Tendo em vista se tratar de uma máquina síncrona trifásica, notam-se os enrolamentos do estator, 
em número de três, o enrolamento de campo e os enrolamentos amortecedores. 

Em todos estes enrolamentos aparecem os efeitos da indução própria e da indução mútua, uma vez 
que eles se interagem simultaneamente. 

De uma maneira geral, podem-se equacionar os efeitos dos fluxos enlaçados segundo a relação 


matricial seguinte. 


ha | |laa lab lacilar LaD taQ || ia 
4 | |tba Lbb Lhe; LbF LbD LbQ || Íb (015) 
to | Jica Lob Lee; LoF LeD LeQ ||io 
AF | |lra Lrb Leci trF LD LrQ|iF 
AD| |LDa LDb LDciLDF LDD LDQ|iD 
12] [Loa Lob tlocitor top toallia 


Onde: 
Aab, Da: flUXOS instantâneos enlaçados nos enrolamentos indicados pelos índices 
Lab,F Do: indutâncias de estator e rotor por fase 


lab, FD,Q: Correntes instantâneas de estator e rotor por fase 


Considerar para a designação da indutância L; que: 
j 
j=i: para a indutância própria do enrolamento considerado; 


jai: para a indutância no enramamento i devido ao enrolamento j. 
A equação matricial 015 anterior pode ser entendida como a representação linear entre correntes e 
fluxos o que confere linearidade nesse processo de equacionamento da máquina síncrona através 


de suas indutâncias. 


2.4. Indutâncias Próprias e Mútuas do Estator 

Para obtenção das fórmulas das indutâncias próprias e mútuas do estator de uma máquina síncrona 
considerar os cálculos cujos desenvolvimentos podem ser vistos na referência 6. 

Simplesmente, as expressões são as seguintes: 


Laa = Ls + Lm cos 20 (016) 


Onde: 


1 


FR e | 
E) a 2 
e 

2 SUCH 
Lm=N 
m a 2 


Nas equações anteriores considerar: 

Na : número de espiras da fase a do estator; 
3a: permeância magnética do eixo d do rotor; 
3q: permeância magnética do eixo q do rotor. 


Portanto, Laa possui uma parcela síncrona Ls e uma parcela dependente da velocidade dupla .,, cosze 


Considerando as outras duas fases do estator da máquina síncrona, têm-se as seguintes expressões 
para as indutâncias próprias do estator: 


Lbb = Ls + Lmcose(o-*) (017) 
Emos toada cosa(6 +) (018) 

Para a indutância mútua tem-se que: 
Lab=Lba=-Ms-Lmcos2(9+5) (019) 


Onde: 
Ms = NÊ Jd +3g 
4 
e 
2rSd=sg 
EmenNge = 
m a 2 


Portanto Lap ou Lba possui uma parcela síncrona Ms e uma parcela dependente da velocidade dupla 


Lmcos2(6+5)- 


Utilizando da analogia com a equação 019 relativa à influência da fase a na fase b, considerando-se 


que as outras duas fases do estator da máquina síncrona devem receber tratamento matemático 
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semelhante, já que o sistema é considerado trifásico e equilibrado, tem-se para estas fases as 


seguintes expressões para as indutâncias mútuas do estator: 


Lhe = Lob = Ms -Lmcos2(6-*) 
Lca = Lac =-Ms -Lmcos2(6+*E) 


2.5. Indutâncias Próprias e Mútuas do Rotor 


(020) 


(021) 


As indutâncias próprias e mútuas do rotor independerão do ângulo 6 de giro da máquina síncrona, 


referenciado na equação 002, uma vez que seus valores já se encontram vinculados aos eixos de q 


relativos ao sistema de eixos já em movimento na velocidade síncrona. 


Assim, para as indutâncias próprias do rotor, tem-se, considerando no rotor os enrolamentos 


amortecedores de eixo direto e de eixo em quadratura D e Q, respectivamente, e o enrolamento de 


campo apenas no eixo direto, enrolamento F, seguindo a designação mais comum na literatura. 


Para o enrolamento de campo tem-se que: 


Lrr =Lp 
Para o enrolamento amortecedor de eixo direto: 

LDD=LD 
Para o enrolamento amortecedor de eixo em quadratura: 

Lag =LQ 
Para as indutâncias mútuas tem-se que: 
para campo e o enrolamento amortecedor D: 

ED DE TME 
para campo e o enrolamento amortecedor OQ: 

Lrq=tgr =0 
para os enrolamentos amortecedores De Q: 

LDQ=LtoD=0 


(022) 


(023) 


(024) 


(025) 


(026) 


(027) 
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Nota: Nas equações 026 e 027, as indutâncias mútuas entre os enrolamentos De Qe entre os 
enrolamentos F e Q, nesse sentido, são nulas, isso porque esses enrolamentos fictícios perfazem 


entre si um ângulo teórico de 90º. O fluxo indutor da máquina síncrona é dado pela fórmula 


fo) -[ B.hds, Uma grandeza escalar do cálculo vetorial onde o integrando é o produto escalar do vetor 
s 


': Ea EE ob , das 
indução magnética g pelo vetor relativo à área ds =nds, um vetor normal à superfície S por onde o 


fluxo passa. Essa integral, para o presente caso, é nula uma vez que o ângulo envolvido vale 90º. 


2.6. Indutâncias Mútuas Entre Estator e Rotor 


Considerando que a força eletromotriz, fem, gerada no enrolamento de campo apresenta apenas a 
componente de eixo direto, torna-se possível a obtenção das expressões para as indutâncias mútuas 


do estator com o enrolamento de campo. 


Lar = Ne23d cosg (028) 
ou simplesmente 
Lar = Lra = Mpcos6 (029) 
Onde 
ME=NE 3a 
Analogamente, para as demais fases, tem-se que: 
Lbr = Lrb = Mpcos(0- au ) (030) 
Lc =LFc= Mpcos(os E ) (031 ) 


Para as indutâncias mútuas entre as fases do estator e o enrolamento amortecedor de eixo direto, 


tem-se, considerando a fase ae o enrolamento D: 
LaD = LDa = MD cos6 (032) 


Onde: 


Mp = No? 3d 
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Analogamente para as indutâncias mútuas entre as demais fases do estator e o enrolamento 


amortecedor de eixo direto, tem-se: 


LbD=LDb=MD cos(a-) (033) 


LcD = LDc = MD cos(0+ o (034) 


Para as indutâncias mútuas entre as fases do estator e o enrolamento amortecedor de eixo em 
quadratura, Q, tem-se um caso análogo ao anterior, porém com a projeção executada sobre o eixo 
q. 


Assim 
Lao = La = Mo sine (035) 
Lbo = Lob = MO sin(o- E ) (036) 
Lea = Loc = MQ sin(o+ gu ) (037) 
Onde: 
Mo = N3q 


2.7. Equações dos Fluxos Enlaçados e Desenvolvimento 


Considerando-se a equação 015, relativa aos fluxos enlaçados envolvidos com a máquina síncrona, 


e em seus valores instantâneos, numa forma simbólica mais compacta, vem: 
pesa is fas 
4F,DQ] [LRE LRR] LFD,Q (038) 


Multiplicando os membros da equação 038 pela matriz seguinte indicada na equação matricial obtém- 
se que: 


oca oiee! iso 
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Onde: 


| é uma matriz unitária na dimensão 3x3 e P a matriz da transformação d-g-0. 


Promovendo as operações, inclusive a multiplicação do segundo membro da equação 040 pelo 
produto unitário matricial P'P=1, vem que: 


Elo A E E plo H lisa 
*rDaQ] LO illlre trajL o JJ MlrDal (040) 
E também: 

pi E o e es 

4rDQ| LO lltpe LarlL o J|lFDO (041) 


A equação matricial 041 encontra-se numa forma compactada. Descompactando-a obtém-se a 
seguinte equação: 


Allo 0 0:00 o og 
Ag a? 0 0 Lg 0 0 KMg ig 
AF O kMp O: Le Mp Ollir 
Ap (9) KMp 0 Mp Lp 0 'D 
o) 0 0 kMg 0 0 Lg 


2.8. Equações das Tensões 


O esquema do diagrama da figura 3 mostra os circuitos dos enrolamentos componentes de uma 
máquina síncrona real operando como gerador. 

Considerando o esquema da figura 3, segue então a equação matricial compactada das tensões 
envolvidas na máquina síncrona, em seus valores instantâneos. O ponto sobre a grandeza indica 


derivada desta em relação ao tempo. 
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v=vp+ri+A=0 (043) 


De maneira expandida, a equação matricial anterior toma a seguinte forma: 


| acto a 
va ra o 0:00 0 Ollia ' 
vb o mo;jo o ollip| |M (044) 
da [E ndo (e A AR RR DR PRE 
“vp Loo 0 dm 0 Oi 
0) o 0 0:0 mp Ollipl | 
0 o 0 0:0 0 rallia| |1D 
L no 
'F 


Figura 3 — Enrolamentos Reais da Máquina Síncrona 


Na figura 3 considerar, além das variáveis e parâmetros conhecidos e já definidos anteriormente, a 
seguinte legenda: 

Faabboc FD: resistências dos enrolamentos; 

vr : tensão instantânea de campo; 

in: corrente instantânea de neutro 

In: resistência de neutro do estator; 


L, : indutância de neutro do estator. 
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Para determinação da matriz vrelativa às correntes envolvidas tem-se: 


Í 
rn rn rnllia Ln Ln Ln : (045) 
vnS|rn rn enllibl |Ln Ln Lnlib 
rn rn enllic] |lLn ln ln 


De maneira mais resumida, a equação 045 toma a seguinte forma mais compacta: 


Vn=-Rniab,c—Lniab,c (046) 


De maneira compacta, a equação 043 poderá ser reescrita da seguinte forma: 


E [ta |i[Pade 0 | e | e -0 (047) 


vVEDQ| LO O ReDQI|IFDOQ||* 
AF,D,Q 
Aplicando a transformação d-g-0 à equação 047 vem que: 
E j Vabe | | Yoda (048) 
0 1YFDQ]| |YFDOQ 


Considerando, a partir daqui, toda a dedução algébrica que pode ser vista nas referências 1 e 6, 
obtém-se a seguinte equação matricial das tensões instantâneas em função das correntes 


instantâneas relativas aos enrolamentos da máquina síncrona. 


[volt 0 010 0 0 Tio] 
vd 0 mp wlg: O 0 wkMo ig 
Ru PDS Ra O Pd lo RR A 
2VE Õ 0 CO tr 0 0 “ir 
vp=0 0 0 0: 0 rD 0 ip 
vo=0| | 0 A 0 ro ig 


Lo+3Ln O 0:o 0 o To 
0 Ld O :KME KMp O id 
0 0 lg. 0 0 KM 
0 KMp (0) Lr Mp 0 o 
0 KMp O Mp Lp Ol (049) 
0 o KMg: 0 0 tg |!D 
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Utilizando a equação 049, pode-se representar o circuito equivalente para uma máquina síncrona, 
como se vê na figura 4, com ra= rp= re=r, uma vez que todo projetista de máquina síncrona persegue 
esta condição. 

Pela figura 4, pode-se perceber, no circuito dos eixos d e q, a presença de duas tensões que são 
proporcionais ao movimento atual de rotação do conjunto e que são chamadas tensões de rotação. 
Portanto, as relações obtidas não são lineares, mas, por aproximação, tomando-se a velocidade de 
rotação w como constante, uma vez que a máquina síncrona dispõe um regulador de velocidade, o 
problema fica resolvido obtendo-se a linearização procurada. Lembrar que é mais prático e 


conveniente se trabalhar com equações lineares. 


"e + 3h 
+ . À 
El) Eno L0+3Lnéllo fi 
Eq 
KMF 
'b MR r 
e. KM Ld id Va 
- + 
Ag 
"q Fr 


iq Lq << K Mo —= Lg 


Figura 4 — Circuito Equivalente da Máquina Síncrona 
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2.9. Valores em Por Unidade no Estudo da Máquina Síncrona 


Cabe aqui uma observação, para a utilização das equações da máquina síncrona em por unidade há 
necessidade de se proceder à normalização das grandezas envolvidas no processo. 

Para a normalização das grandezas utilizadas deverá ser adotado um referencial baseado no rotor, 
proposição feita por alguns autores apesar da escolha depender de certos critérios. 

Os fluxos enlaçados mútuos dependem da parcela mútua das indutâncias que compõem o sistema 
magnético da máquina síncrona. 

Considerando as parcelas do fluxo mútuo e do fluxo de dispersão, por definição, pode-se escrever a 
seguinte expressão para o valor total da indutância de eixo direto da máquina síncrona, primeiro caso 


a ser visto e a ser seguido daqui para frente neste tópico. 


Lg=Lma+Ldd (050) 
onde 


Lma : indutância mútua útil referente ao eixo d; 


Lmo : indutância de dispersão útil referente ao eixo d. 


A indutância de dispersão tem a ver com o fluxo perdido no entreferro. 
Assim, ficam definidas, para todas as indutâncias em jogo, as relações seguintes referentes às 


parcelas úteis de cada uma em seu respectivo enrolamento conforme é mostrado na figura 4: 


Lma =La Lad 
Lmg= Lg Log 

Lmp=LF-LdF 
LmD=LD-LdD 


Imq=ta-Ldo (051) 


As relações elementares para um determinado valor de fluxo enlaçado As, aqui tomado como um 
valor base que construirá a base de valores a ser escolhida no processo de normalização envolvendo 
outras grandezas tais como a tensão base Vs, a indutância base Ls, a corrente base |; e rotação 


base ws, são as seguintes: 


(052) 


48-=Lp!B 


(053) 


vB 
AB=—2=LBl 
B B'B 
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com 


1 (054) 


Notar que a equação 054 define uma grandeza que não deve ser confundida, trata-se do tempo base 
que é exatamente o inverso da velocidade angular elétrica, o que para o sistema elétrico nacional, 
SIN, tem seu valor aproximado para wp= 2760 = 377 rad's . 


As seguintes premissas serão obedecidas: 
sp=VpIp > potência para sistemas monofásicos; 


sp =3Vp1p potência total para sistemas trifásicos. 


Onde: 

Se : MVA monofásico de base para sistemas monofásicos; 

Se: MVA trifásico de base para sistemas trifásicos. 

Podem-se então definir, a partir de certas condições, as grandezas básicas do rotor. 


Assim, tem-se que: 


Amd = Lmd!B = KMF IrB = KMD IDB 


AmF = KMF IB = LmF lFB = MRIDB 


AmD = KMp IB = MR IFB = LmDIDB 


AmQ = KMg !B = Lma!aB (055) 
Nas equações 055, a idéia foi obter uma escolha de correntes base que, a partir da corrente base 
circulante por um determinado enrolamento, através de sua indutância própria útil, fosse possível 
obter uma corrente induzida base nos outros enrolamentos por esse primeiro enrolamento 
influenciado por intermédio das indutâncias mútuas. 
Exemplo 1 


Utilizando-se das equações 055, provar que as indutâncias mútuas em pu KMFue KMpu são iguais. 


Solução: 
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Das equações 055, podem-se escrever as seguintes relações: 


Amd = Lmd !B = KMF IFB = KMD IDB 


Tomando-se a última igualdade tem-se: 


KM, ) 
KM IFB = KMpDIDB -. E = E (056) 
Porém, sabe-se ainda da equação 053 que: 

meg =—VB ea evÊ (057) 

wBIFB wBkMeB 
vB vB (058) 

KMpB = RA Inpp=——P 

DB wBIDB DB wBkKMpB 


Levando na equação 056 as expressões encontradas nas equações 057 e 058 para as correntes de 
base Is e Ipe, obtém-se: 


V 
B 
wBkMpB 


Por analogia ao resultado obtido anteriormente, pode-se generalizar para os circuitos dos eixos d e 


q da máquina síncrona que: 


Lmau = tmFu= tmDu= Meu = KMpy = May 


Lmqu = tmQu=KMou (059) 


As equações 059 são de extrema importância quando da elaboração dos circuitos equivalentes de 


eixos de q, em pu, da máquina síncrona. 


Nota a respeito do exemplo 1: Observar sempre que, na definição dos valores base para as 
indutâncias KMF e KMpe, nas equações 057 e 058, a tensão base se refere ao enrolamento em 
questão, porém a corrente base se refere ao enrolamento que criou o fluxo magnético. Essa 
proposta de base de valores é assim adotada para que os valores em pu das indutâncias mútuas 
sejam iguais, uma particularidade opcional visando à obtenção de facilidades para futuros cálculos. 
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Exemplo 2 


Comprovar as equações 055 a partir de uma condição genérica. 


Solução: 

Para a comprovação propriamente dita, será utilizado nas variáveis um subscrito, apenas nesse 
exemplo, para que fique bem indicado a qual dos circuitos equivalentes d, D, F, q, Q o cálculo se 
refere. 


imaginando uma escolha genérica para as correntes de base circulantes, tem-se o seguinte 


raciocínio. 


a) enrolamento d 


Os fluxos enlaçados vistos do enrolamento d serão os seguintes então: 


Amdg = tmd!B 
(060) 


AmFg = LmFIFB 


AmDg = tmD!DB 


Para os três fluxos anteriores, a tensão base será sempre a tensão referente ao eixo d, ou seja, Vs, 
já que o efeito se dá no enrolamento de estator. 
Tomando-se os valores em pu para cada um dos três fluxos enlaçados da equação 060 vem que: 


Amdy - LmdlB — Lmd (061) 
AmdgB ImaB!B LmdB 


Amdau = 
A indutância base referente ao eixo d é dada por: 


va (062) 


Com a expressão anterior, a equação 061 pode ser assim reescrita: 
(063) 


wB! 
Amdgu ar Lmd 


Por analogia, os fluxos enlaçados para Amrdu € Ambdu Serão: 
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vBIDE kmp (065) 
4: 


Para que, por suposição, os fluxos enlaçados em pu, relativos ao eixo em questão, sejam iguais, as 


seguintes relações deverão ser cumpridas: 


Amdgu = AmFgu = AmDgu (066) 
Deste modo, tem-se então que: 
IB Lmd = !FB KMF = DB KMp = Amd (067) 
b) enrolamento q 
Estendendo o raciocínio para o enrolamento g vem: 
IB Lmg = !QB KMQ = Amg (068) 


c) enrolamento F 


Para o enrolamento de campo, em termos de fluxo magnético enlaçado, considerando-se a 


circulação da corrente base desse enrolamento, tem-se que: 


AmFE = LmFIFB 


Tomando o valor em pu para fluxo magnético enlaçado desse enrolamento, obtém-se: 


— LmFIFB (069) 


A 
DÉFU  LmFBlFs 


A indutância base referente ao eixo F é dada por: 
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VrB (070) 


LmFB= REIS 
Com a expressão anterior, a equação 069 pode ser assim reescrita: 
inf SEI Lp (071) 
Por analogia, os fluxos enlaçados para Amcru € Ambdu serão: 
mdFy = es KMF (072) 
e 
AmDey = vBID MR (073) 


VrB 


Também aqui, para que, por suposição, os fluxos enlaçados em pu, relativos ao eixo em questão, 


sejam iguais, as seguintes relações deverão ser cumpridas: 


AmdFy = AmFey = AmDFy (074) 


Deste modo também, tem-se então que: 


IrB LmF = !B KMF = IDB MR = AmF (075) 


d) enrolamentos De Q 


Por analogia, podem-se escrever as seguintes equações para os fluxos enlaçados nesses 


enrolamentos: 


IDB LmD = !BKMD = IFB MR = AmD (076) 


IB LmQ = |B KMQ = AmQ (077) 
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Ficando, deste modo, comprovadas as equações 065. 


Nota a respeito do exemplo 2: Multiplicando-se as equações 067, 075 e 076, relativas ao eixo d, 


respectivamente, por /s, Ir e los, vem que: 


Amd!B =Lmd!B? =KMFp IrBlB=KMpIDB!B (078) 

e 
AmF FB =KMF IB lFB= LmFlFB2=Mp !FBIDB (079) 

e 
AmD!DB = KMp!B Ipe = MR !FB!DB = LmDIDB? (080) 


Do mesmo modo para o eixo q, multiplicando-se as equações 068 e 077, respectivamente, por lp e 


los, vem que: 


Amg!B=Lmg 182 =KMa las !B (081 ) 


AmQ!QB =KMa IB las = Lmo los? (082) 


Por uma imposição do problema, as equações 078, 079 e 080, relativas ao eixo d, são iguais, logo 
vale a seguinte relação primordial: 


LmalB? =LmFlFE? = LmD!DE? = KMF IFB!B = KMp IpB!B = MR !FBIDB (083) 


Do mesmo modo para o eixo q, as equações 081 e 082 são iguais, logo vale a outra relação 


primordial: 
Lmg!Bê =Lmo!aBÊ = KMa los IB (084) 


As equações 078, 079, 080, 081 e 082, chamadas aqui de equações primordiais, são relações que 
confirmam que as potências desenvolvidas nos diferentes enrolamentos dos eixos de q da máquina 
síncrona, em termos dos valores base escolhidos nesse contexto, são iguais. 

Na realidade, as equações 055 são uma consequência da condição imposta relativa à conservação 


da energia entre os enrolamentos dos eixos de q da máquina síncrona. 
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Considerando as equações 083 e 084, além de se manter a mesma potência base de valores, podem- 


se escrever, para os três enrolamentos do rotor, as seguintes expressões: 


re Ve [ma ima MME KMp j (085) 
B VrB Vime KMF Imp Mp 


DB VB 2 Emo ctmo AMD KMF cs (086) 
IB VDB VimD KMD LmD Mp 


IoB VB Emg ima  KMQ Ro (087) 
É 


B Vos Vime KMo Lma 


Assim, as variáveis de corrente de base e de tensão de base são escolhidas de modo que em 


qualquer enrolamento a potência base seja sempre a mesma. 
Portanto, conhecendo-se Ls e Rs, das bases do estator, vem: 


RrB. 1 RDB. 1 Ros 1 (088) 
RB Kp2 “RB Kp2 “RB Ko? 


LB. 1 LDB. 1 los 1 (089) 
LB KF2 “LB Kp2 is Ko? 


E para as indutâncias mútuas será utilizada uma propriedade extraída da equação 085, senão veja- 
se. 


(090) 


Lmd KMF 2 
= - KMpº=LmgL 
KME LmF md -mF 


Nota: Nessa propriedade, nota-se que, considerando o efeito mútuo entre dois circuitos 
magneticamente conectados, o produto das indutâncias próprias úteis de cada circuito elétrico é igual 
ao quadrado da indutância mútua verificada entre ambos. Deste modo, por consequência, levando- 
se em conta uma base de valores previamente definida, o quadrado do valor da indutância mútua 
base entre dois circuitos magneticamente conectados será igual ao produto dos valores das 
indutâncias próprias base de cada um dos circuitos envolvidos. 


Assim: 
(091) 


MrB . LFB 
LB Mpb 


MrB2 =LBLFB -. 
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Mas, da equação 085, tem-se que 1, Kc eM valores atuais o que permite escrever em valores base 
Lmd - kr 
KMF 


que: 


ME ad (092) 
LB KF 


Mas, da equação 086, tem-se que Lmd - Kp EM valores atuais o que permite escrever em valores 
KMp 


base que: 


MDB 1 (093) 


Mas, da equação 087, tem-se que Lmg Ko EM valores atuais o que permite escrever em valores base 
Pads 1 Q 
KMa 


que: 


MoB 1 


' (094) 
LB Ka 


Da equação 085 tem-se que: 


KMp ur imakMD po. MA 1 KMp (095) 


MR — Lmd MR “Lmd KF Lma 
Da equação 086 tem-se que: 
kMp 1 (096) 
Emd KD 


Logo, a partir da equação 095, obtém-se para a equação 096 que: 


MR 1 (097) 
Lmd KFKD 
Em termos de relação base, a partir da equação 097, vem que: 
MpB. (098) 
LB “KFKD 
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A equação 098 completa o conjunto de equações envolvendo os valores base fornecendo o caminho 


para o cálculo das demais relações. 


2.10. Equações em Por Unidade 


Colocando a equação matricial 049, que se encontra em valores atuais, em por unidade obtém-se 


que: 
VB Vou o+8m O 0: 0 0 0 Il iBiou | 
VB Vau 0 r ug; 0 0 wkMg|| IBidu 
Vevgu || O culg  r icwkME -wkMp o O | IBigu || 
-VFEBVFu 0 0 O LF 0 0 IFBIFu 
0 0 0 0: 0 [D) 0 IDB'Du 
0 0 Do pop .dO 0 ro |loBiou 
: “| IBiou 
Esta so GR vês. . 
IBidu 


0 Lg O iKME KMp O : 
RO qui ERG DO ao MO EQ 

0 kMp 0 : lr Mp 0) 

0 KMp 0 Mp Lp 0 E 

0 0 KMg: 0 O tg |lpBiDu 


(099) 


Analisando cada uma das tensões do vetor tensão, em separado, iniciando-se pelo eixo d, vem que: 


IB. IB. IQB . Bº 
Vdu tr iduy + wWLg — ioy + WKMo = iQu + Lg —— idu+ 
du Va du q VB qu Q vB Qu ora du (1 00) 


e 
+ KM 'EB ;ry+ KMp'DB ipy =0 
vB vB 


Sabe-se que + 1, assim a equação 100 se transforma na equação 101. 
LeWp Ve 


Vdu + ruidu + Wulquigu + wuKMQuiQu = 
(101) 


iDu=0 


Lay * KMpyº KM, 
AU ;nuA EU j py Du 
WB WB WB 


Analogamente, para o eixo q, vem: 
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Para o circuito de sequência zero vem que: 


Vou +(r+3m É ou +(Lo+3Ln) 
B 


Assim, obtém-se a seguinte expressão: 


r+3r, Lo + SL 
onde aelttn O pot 
vB vB 


Para o circuito de campo vem: 


E, para os enrolamentos amortecedores, tem-se: 


Vqu — Wulduidu + tuigu — VukKMFuiFu — vukKMDuiDu + 


IBwB 
VBwB 


lou =0 


; Paris 
vou + Rou iou + -2É iou=0 
WB 


iFu- 


: KMpy * L M E 
VFu-fFuiFu Er Idu E E iDu=0 


e e e 
idu+ Mau iFu+ LDu iDu=0 
WB WB 


KM 
rDu iDu + ——Du 
WB 


; KMqu *º a 
"Qu'Qu + E igui e iQu=0 


(102) 


(103) 


(104) 


(105) 


(106) 


Colocando as equações 101, 102, 103, 104, 105 e 106 na forma matricial completa, e normalizando 


o tempo nas parcelas que apresentam derivadas, obtém-se que: 
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Lou O 010 0 o lo 
0 Ly O :KMpy KMpy O idu 
1|0 0 tLuio 0 KMau| ligu 
) 
0 
0 


wB KMpy O ley Mau 0 E 
KMpy O iMpu lou O ||ifu 
0 KMqu: O 0 tou ||iDu (107) 


2.11. Circuito Equivalente da Máquina Síncrona 


Como condição a ser obedecida nos desenvolvimentos que se seguem, deve-se levar em conta que 
o sistema é equilibrado eletricamente com io = O. A base de valores escolhida é aquela descrita no 
item 2.10, precedente, cuja escolha partiu de certas premissas a fim de se fazer com que as 
indutâncias mútuas fossem iguais em seus valores em pu. 

Os fluxos enlaçados da máquina síncrona, assim como foi visto no item 2.9, podem ser relacionados, 
em pu, com os fluxos mútuos e os fluxos de dispersão, respectivamente, de modo que se pode 


escrever de maneira genérica que: 
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ATu = Adu + Amu (108) 


ou seja, o fluxo enlaçado total Ar, de um determinado circuito magnético, é a soma dos fluxos de 
dispersão Aq e mútuo Am, desse mesmo circuito. 


Assim, tem-se que: 


ATdu = Addu + Amdu = Ldduidu + Lmduidu + KMF uiFu + KMDuiDu 


ATqu = Adqu + Amqu = Laquigu + Lmquigu + KMQuiQu 

(109) 
ATFu = AdFu + AmFu = Laruiru + EmFuiFu + KMpuidu + MRuiDu 
ATDu = AdDu + AmDU = LaDuiDu + LmDuiDu + MRuiFu + KMDuidu 


ATQu = AdQu + AmQu = LaQuiQu + LmQuiQu + KMQuigu 


Das equações 059, sabe-se que: 


Lmdu = LmFu = LmDu = KMpu = KMDU = MRU 


Lmqu = tmQu = KMqu 


As duas equações anteriores permitem que se escrevam as seguintes relações em pu: 


Lmdu = KMFy = KMpu = Mu = Lidu (110) 


Lmqu = tmQu = KMQu = Ligu (111) 


E criando os fluxos enlaçados associados a cada uma das indutâncias gerais das equações 110 e 
111 vem que: 


Aidu = Lidu (idu +iFu +iDu) (112) 


Aigu = Ligu(igu + iQu) (113) 


O significado físico das equações 112 e 113 pode ser representado considerando-se cada um dos 
eixos de q como está mostrado na figura 5. 


Para o circuito equivalente completo, devem-se levar em conta as equações em tensão 107, já 


considerada a condição de se estar trabalhando com um sistema equilibrado trifásico com as 
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resistências das fases do estator tomadas com o mesmo valor, assim como se tem procedido até 


aqui. 


Para o eixo d pode-se escrever que: 


; 1º 
Vdu ! Fu ldu + wB Adu+ WuAgu =0 (1 14) 


Da primeira equação das equações 109 extrai-se a expressão para Aw total, o que leva a seguinte 


expressão: 


e * ae » Y e 
ATau = Leiau idu+ Lmdu idu+ KMpu iFu+KMpu iDu (115) 
Assim, a expressão da tensão 114 para o eixo d tomará a seguinte forma: 


116) 


; 1 e Lia Ma, 45 
Vdu + tu idu + — Ladu i du+ WuAgu + idu (idu+ iFut Du) o( 
wB WB 


De maneira similar, chega-se às equações de tensão para os demais enrolamentos. 


Para o enrolamento de eixo q tem-se que: 


; 1 Liqu,.* * 117 
Vou + ty igu + —L i WuAdu + iqut =0 ( ) 
qu + "ul!gu wB dqu 'gu” YuAdu op ( qu Qu) 


Para o enrolamento de campo tem-se que: 
(118) 


; 1 - Lidu, 0.88 
=Vpu+rFulFut——LdFu i Fut (idut Fut iDu)=0 
WB WB 


Para os enrolamentos amortecedores De Qtêm-se: 


: 7 me L; RA Ea a 
vDu =0:: DuiDu + LaDuiDu op (idu” iru+ipy)=-0 (119) 


1 * Hiqu 
VQu = O. "Qu iQu *og dou iQui DE (igul iQu)=0 (120) 
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As figuras 6 e 7 mostram os circuitos equivalentes de eixo d e eixo q, completos em tensão, 


respectivamente. 


Figura 7 — Circuito PU Equivalente de Tensão da Máquina Síncrona para o Eixo q 


2.12. Autoexcitação em Máquinas Síncronas 


O fenômeno da autoexcitação em máquinas síncronas é um processo que pode ocorrer quando a 
máquina estiver conectada a um circuito elétrico de característica capacitiva. 

Deste modo, qualquer tipo de dispositivo capacitivo que esteja acoplado à máquina síncrona pode 
vir a provocar esse fenômeno. 


Para análise, será suposto o sistema mostrado na figura 8. 
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Máquina Carga 
Capacitiva 


TESS C 


Figura 8 — Circuito Máquina Síncrona Conectada a uma Carga Capacitiva 


Figura 9 — Diagrama Vetorial para o Circuito Máquina Síncrona Conectada a uma Carga Capacitiva 


No diagrama da figura 9 podem-se observar as grandezas envolvidas no processo sendo ia corrente 
circulante, quando o circuito aposentava apenas uma determinada carga indutiva, e sendo ia 
corrente capacitiva adiantada da tensão terminal v após a entrada da carga capacitiva. 


O seguinte equacionamento básico é possível. 


va =Vsin(-89)=—v sind 
vq =vcos(-89) =vcosdg 
Com 
ig =icsin(m/2-00)=i, cosdo -. cosdg =iy a 


ig =ip cos(m/2-00)=iç sinôg -. sindg =ig fic 
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Considerando-se apenas a presença da carga capacitiva vem que: 


v=-j Xcic or v=Xeiç 


Nota: a grandeza tensão v foi representada anteriormente na forma vetorial para indicar sua 


posição nos eixos d-g. Assim, as tensões va e vytomam a seguinte forma: 


q Vis ; : 
Vg Na ta Kolig=-Xcig 


(121) 


ig Vo. : ; 
va VC L)=(D)g)=Ko)ig =Xcig 
te te 


Onde: paso esci=) 
wC 


Considerando o fluxo enlaçado pelo eixo d, através das equações 109, desconsiderando o 
enrolamento amortecedor, tem-se que: 


Adu = Lduidu + KMFuiFu (122) 


Para o fluxo enlaçado pelo campo, através das equações 109, desconsiderando o enrolamento 
amortecedor, tem-se que: 


Aru = LruiFu + KMpuidy (128) 


Multiplicando a equação 123, membro a membro, pela relação «mr, vem que: 
Lru 


2 
Ma SRA MB q (124) 
Lru Lru 


Considerando-se as referências 1 e 6, as seguintes relações são definidas: 


Eu nr = 5Equ = equ que e a fem instantânea da maquina sincrona para a condição operaçao em 
Fu 


vazio; 


KM | 
EU vpy = KMryiru = 3 Eta = efa 


"Fu 


que é a tensão instantânea de campo da máquina síncrona; 


x (KMpu 2 que é uma relação primordial. 


Lgu-L'du E 
7 
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A partir das relações primordiais anteriores, a equação 124 tomará o seguinte aspecto: 
egu = KMFuiFu + Ldu idu — L'du idu (125) 
Considerando-se agora a equação 123, a equação 125 poderá ser reescrita da seguinte forma: 
equ = Adu- L du idu 
Passando a equação anterior para valores atuais, e considerando que w= we, vem que: 
eq =wAg-wL'gid (126) 


Das equações 101 e 102 obtêm-se as seguintes expressões para as tensões va e vy, desprezando- 


se as resistências envolvidas. 


Fa (127) 
va = Ag 
Das equações 121, a partir das equações 127, conclui-se que: 
wAg =Xcid (128) 
wAg =Xcig 
Assim, a equação 126 tomará a seguinte forma: 
eg =Xcig-wLgig (129) 
O que, finalmente, conduz ao seguinte valor da corrente ig: 
j a E (130) 


Se o mesmo raciocínio for levado em consideração para um determinado enrolamento de eixo q, 
onde todas as condições forem respeitadas, a seguinte expressão para a corrente iq será obtida: 
UC (131) 
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A equação 104 fornece a expressão em por unidade da tensão de campo, onde a resistência de 
campo é considerada, o que permite que se escreva a seguinte expressão, considerando-se o fluxo 


através do campo da máquina síncrona: 


vrufFuipu= AFu =0 (132) 


Na equação 132, mais precisamente no termo da derivada do fluxo, a parcela 1/ws foi incorporada 


à variável Ary, O que garante estar a variável tempo sendo dada em segundos. 


O valor do fluxo enlaçado no campo da máquina síncrona será dado por: 


t 
Aru = "Fu od = iruddt (133) 
u 


Multiplicando a equação 133, membro a membro, pela relação «me, vem que: 
Leu 


KM, reu ft v 
EU ypy = EU [ (kMpy-ÉU  kMpyipudt 
Lru Lry"O fFu 


Considerando as relações primordiais e sabendo-se da referência 1 que rião= “Eu é a constante de 
Pina a 


"Fu 


tempo transitória de eixo direto a circuito aberto do estator da máquina síncrona , vem que: 


Ea is 134 
egu = Tao lo!Stau — KMpy iry)dt (134) 


kMFy iFu» da equação 125, obtém-se a seguinte 
Trazendo na equação anterior 134 a expressão para 


expressão mais geral: 


ad no (135) 
equ Tao lo! etau [equ — (Ldu — Lu Jidu] Jdt 
Passando a equação 135 para valores atuais obtém-se que: 
pet La- L'g Vig] dt 
eg Ta loteta leg-w(Lg-L'g dig]! 
Ou ainda: 
e ps ago Vig] (136) 
w Tdaodw uw 
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Trazendo na equação 136 a expressão da corrente ia, da equação 130, tem-se que: 


q 1 ste e , wC 
lot E gt —S— egdt 
w Tão tw u t=w LG 


Ordenando a equação anterior segundo cada membro e passando o resultado para o plano complexo 
da transformada de Laplace, sendo S operador complexo, obtém-se que: 


1 a 1 wC 
ea(s) ==—=[—efd(s)— —eg(s) + (Lg —L'd ——s — eg(S] 
US) RS (9528) RA 


Rearranjando-se a equação anterior vem: 


14 Lg =L'g JC 

eusli+ (Lg-L'd) . 

q TãoS T 2p TãoS 
w wTdoS TywSi-wLgo Tão 


erd(s) 


Multiplicando a equação anterior pelo fator wrmws obtém-se que: 


, ; w2C 
eg(s)N+TdoS-(La-L'd -—s J=efd(s) 
1-wíL'gC 


As operações subsequentes permitem que se chegue à função de transferência “a(s). Assim, 


efd(S) 
resolvendo a equação anterior, tem-se que: 
1-w2LgC 1 
a(s) 27 nl 7 fds) + a(o) 
Tão! -w L'gC) “do (137) 
Para condições inicias nulas, vem que: 
EA 
Pa(S) Tdo 
HAS os 1-w2LgC (138) 


Tão (!-w2L'4C) 
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Se for considerada a existência de condições inicias diferentes de zero, eg(0), então a seguinte 
equação representará o problema em questão: 


A equação 138 traz, finalmente, a função de transferência entre as tensões de campo exe a tensão 
de excitação esquando do carregamento da máquina síncrona por cargas capacitivas, considerando- 
se o eixo q. Trata-se de uma função de transferência de primeira ordem com um pólo real cujo valor 
depende dos parâmetros envolvidos no processo. 

É evidente que a combinação de valores entre os parâmetros pode designar a existência de um pólo 
positivo ou um pólo negativo na função de transferência no plano complexo de Laplace o que terá 
uma influência importante na questão da presença ou não de sobretensão no campo. A questão 
torna-se crítica quando o valor do pólo for negativo o que denota, pela equação 138, que a resposta 
no domínio do tempo do sistema será regida por uma exponencial crescente o que representará, 


fisicamente, a presença de sobretensões elevadíssimas. 


Para o caso do eixo d, considerando-se os cálculos efetuados até aqui, a seguinte equação poderá 
ser obtida por analogia ao cálculo anterior verificado no eixo q, equação 136: 


e 1 st, e ni 
dit rê Lg-Laliglat (139) 


Na equação 139, notar que, em relação à equação 136 anterior, a tensão de campo de eixo q não 


existe. 


A solução da equação 139 seguirá o mesmo procedimento que foi verificado na solução da equação 
136. Assim, tem-se que: 


1-w2LgC 
ed(s)[S+——s =] =8d(0) 


Tgo(!-wtLqo) (140) 
Onde: 
1 
O MEMO (141) 
S+ q 


Tai wtbio 
Nota: A questão da solução das equações 137 e 140 tem a ver com a forma que a solução de um 


circuito elétrico admite, no domínio da transformada complexa de Laplace, quando existem condições 


iniciais no sistema em processamento. 
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2.13. Análise da Autoexcitação em Máquinas Síncronas e Conclusões 


Pelo que foi determinado no item anterior é possível estabelecer as faixas para os limites de valores 
dos parâmetros da máquina síncrona em função da questão da sobretensão originária no fenômeno 
da autoexcitação provocado pela interação entre os parâmetros de máquina versus os parâmetros 
relativos à carga e, por conseguinte, interferindo com a parte de transmissão do sistema. 


Considerando que, para haver sobretensão, o pólo da função de transferência a(s) deve ser 
efd(S) 


negativo, as seguintes hipóteses são consideradas. 
Primeira hipótese, com base na equação 138: 


1-w2LgC<o e i=W Ligo 


I-w2LgC< 0. ar <wLg 
wC 


1=wÉLy OS 0. gi >wL'q 
wC 


Ou seja 
e O 
wLl'g<— <wLg 
we (142) 
Segunda hipótese, com base na equação 138: 
1-wêLgC>0 e t-wêLgC<o 
1-w2LgC>0 E E, 
wC 
I-wêLgC<o : LD cwl'g 
wC 
Ou seja 
1 ; 
wLg <— <wL'g 
we (143) 
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Como se sabe, em termos de máquina síncrona, a impedância de regime permanente ou síncrona 
Xa será sempre maior que a impedância a essa associada de regime transitório X' fato que inviabiliza 


a segunda hipótese. 


Estendendo a análise para o caso em que se foque o eixo d a seguinte condição deverá ser 


verificada: 


1 
wL'g<— <wL 
q wC q 


(144) 


A equação 143 estipula a faixa de valores para a reatância capacitiva relativamente aos valores das 
reatâncias de eixo direto transitória e de regime permanente onde a possibilidade de ocorrência do 
fenômeno de autoexcitação no enlace do fluxo de eixo direto é passível de ocorrer. 

De modo semelhante, a equação 144 estipula a faixa de valores para a reatância capacitiva 
relativamente aos valores das reatâncias de eixo quadratura transitória e de regime permanente para 
O caso em que haveria sobretensão pelo efeito de autoexcitação. 

O fenômeno vem então associado ao crescimento exponencial da tensão de excitação da máquina 
síncrona o que pode trazer danos severos ao conjunto elétrico independente da atuação rápida do 


sistema automático de excitação. 
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3. Considerações Finais 


Pelo que foi exposto, torna-se importante o conhecimento dos parâmetros que vão representar a 
carga que estará conectada à máquina síncrona. Existem situações em que determinados arranjos 
são possíveis em função da abertura de determinadas linhas de transmissão ou mesmo situações 
de baixa carga em que a carga capacitiva torna-se relevante. Nesses casos é interessante analisar 
todos os arranjos possíveis determinando o valor da carga capacitiva equivalente no sentido de se 
evitar o fenômeno da autoexcitação. Outro ponto que também deverá ser observado são os 
esquemas em que são utilizados compensadores capacitivos no sentido de controlar os níveis de 
tensão. Tais compensadores, dependendo do tamanho, podem concorrer diretamente para culminar 


com o fenômeno da autoexcitação. 
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